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Abstract 

We describe the present status of the pion distribution ampUtude (DA) as it originates from several 
sources: (i) a nonperturbative approach based on QCD sum rules with nonlocal condensates, (ii) an 0{as) 
QCD analysis of the CLEO data on F'^^ ^{Q"^) with asymptotic and renormalon models for higher twists, 
and (iii) recent high-precision lattice QCD calculations of the second moment of the pion DA. We show 
predictions for the pion electromagnetic form factor, obtained in analytic QCD perturbation theory, and 
compare it with the JLab data on F-^{Q'^). We also discuss in this context an improved model for nonlocal 
condensates in QCD and show its consequences for the pion DA and the 77* — > tt transition form factor. 
We include a brief analysis of meson-induced massive lepton (muon) Drell-Yan production for the process 
t:~ N — * fi'^ ii~ X, considering both an unpolarized nucleon target and longitudinally polarized protons. 
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INTRODUCTION 



The pion DA parameterizes the matrix element of the nonlocal axial current on the light cone 

1 

(0|J(z)7^75C(z,0)^.(0) |7r(P)) = iUP^ [ dx e"-(^^)(^Tw-2(^^ ^2) 



(1) 



The gauge-invariance of this DA is ensured by the Fock-Schwinger connector [4] (Wilson line) 



C{z,0) =Vexp 







inserted between the two quark fields. The physical meaning of this DA is quite evident: it is 
the amplitude for the transition 7r{P) u{Px) + d{P{l — x)). It is convenient to represent the 
pion DA using an expansion in terms of Gegenbauer polynomials CJ (2x — 1), which are one-loop 
eigenfunctions of the ERBL kernel [s, 4|, i. e., 



n>l 



(2) 



where (f^^{x) = 6x{l — x) is the asymptotic pion DA. This representation means that all scale 
dependence in 937^(2;; /u^) is transformed into the scale dependence of the set {a2(/i^), a4(/u^), . . .}. 
We mention here that the ERBL solution at the 2-loop level is also possible using the same repre- 
sentation 01) 

In order to construct reliable QCD SRs for the pion DA moments, one needs, as has been shown 
in y, [3], to take into account the nonlocality of the QCD vacuum condensates. For an illustration 
of the nonlocal condensate (NLC) model, we use here the minimal Gaussian model 

(g-(0Mz)) = (g-g)e-l-'l^?/8 (3) 

with a single scale parameter = (A;^) that characterizes the average momentum of quarks 
in the QCD vacuum. Its value has been estimated in the QCD SR approach and also on the 

nn 



lattice 
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]J,ll3|,ll4: 



A2 = 0.35 - 0.55 GeV^ . 



(4) 



Let us write down, as an example, the NLC QCD SR for the pion DA ipTj{x). To derive it, one 



starts from a correlator of the currents J^^{x) and Jl^^.f^fiO) = d{0) 7175 (nV)^n(O) with a light-like 



vector n, n-^ = to obtain next SRs for the moments (x^),^, and finally to apply the inverse Mellin 



transform and arrive at {x =^ (Pt^{x). As a result, we then find 





^iix;A) (5) 



i=2V,3L,4Q 




= 0.6 GeV^ 


(b) 
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Fig. 1: Panel (a): The contributions to Eq. ([5]) are shown due to the perturbative loop (dotted line) and the 
four-quark condensate: Lp^^^{x) (standard QCD SRs); (^4Q^(a;, = 0.55-0.80 GeV^) (the NLC QCD SRs). 
Panel (b): Allowed values of the parameters 02 and 04 of the bunches ([6]), evaluated at /i^ — 1.35 GeV^, 
for three values of the nonlocality parameter: = 0.4,0.5, and 0.6 GeV^. Panel (c): The inverse moment 
{x^^)'K-, obtained using the NLC SR ([5|), is shown as a solid line (central value) with error-bars represented 
by dashed lines. 

with A = A^/M^. The local limit A ^ of this SR is specified by the appearance of (5-functions 
concentrated at the end-points x = and x = 1, for example, (^4q(x; A) = 9[5(x) -|- (5(1 — x)]. 

The minimal Gaussian model ([3]) generates the contribution (/74Q(a;;A), shown on the right 
panel of Fig. [1] in comparison with the perturbative one for the standard (local) and the NLC 
types of that SR. We see that due to the completely different behavior of the perturbative and 
condensate terms in the local QCD SR, it is difficult to reach a reasonable consistency. In contrast, 
the NLC contribution behaves similar to the perturbative one. Just for this reason, we have 
a very good stability in the NLC SR case. After processing SR ([5]) for the moments = 
Jq (pn{x) {2x — 1)^ dx, we can restore the pion DA (fTr{x) by demanding that it should reproduce 
the first five moments i = 2 , 4 , . . . , 10, using to this purpose the minimally possible number 

of Gegenbauer harmonics in representation ([2]). It comes out that the NLC SRs for the pion DA 
yield a bunch of self-consistent two-parameter models at fi^ ~ 1.35 GeV^: 

^^LC(^. ^2) ^ ^As(^) ^ ^^(^2) ^3/2(2x - 1) + a4(/.g) '(2x - 1)] . (6) 

The central point corresponds to af^^ = -h0.188, af^^ = -0.130 for = 0.4 GeV^, whereas 
other allowed values of the parameters 02 and 04 are shown on the central panel of Fig. [1] as a 
slanted rectangle [15,]. Because the inverse moments of all these pion DAs equal 

= 3.17 ±0.20, (7) 

being in good agreement with the estimation dictated by an independent SR for this moment, we 
term this bunch self-consistent. This SR can be obtained through the basic SR ([5]) by integrating 
over X and using the weight x'^ (at /Uq ~ 1.35 GeV^) to find 

(x-i)SR = 3.30 ±0.30, (8) 
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Fig. 2: Panel (a): Comparison of the profiles of three pion DAs: BMS (solid line), CZ (dashed line), and 
the asymptotic DA (dotted line). Panel (b): Histograms for the contributions of different bins to the inverse 
moment {x^^).^ shown for the CZ and the BMS DAs. 

shown graphically on the panel (c) of Fig. [TJ 

It is worth emphasizing that the moment could be determined only with NLC SRs 

by virtue of the absence of end-point singularities. Comparing the obtained pion DA with the 



Chernyak-Zhitnitsky (CZ) one 16[], reveals that although both DAs are two-humped they have 
distinct characteristics: the BMS DA is strongly end-point suppressed, as illustrated in Fig. [2^a), 
while the CZ one is end-point dominated. To display this property more explicitly, we show on 
panel (b) of this figure the comparison of the BMS and the CZ contributions to the inverse moment 
(x~^)t^^ grouped in different bins and calculated as ^^^^''^'^ u~'^'^{u) du (normalized to 100%). 



II. ANALYSIS OF THE CLEO DATA ON F^^.^(Q2) 

Many studies Q, 17, 18, 19 



211 ] have been performed in the literature to determine the pion 



22] on the pion-photon transition form factor F^^*^(Q 



DA using the high^recision CLEO data 

In particular, in [.2l| we have used Light-Cone Sum Rules (LCSR) [17|, ll8|] to the next-to-leading- 
order accuracy of QCD perturbation theory to examine the theoretical uncertainties involved in 
the CLEO-data analysis in order to extract more reliably the first two non-trivial Gegenbauer 
coefficients 02 and 04, which parameterize the deviation from the asymptotic expression 99^*^. 

Let us clarify why it is advantageous to use LCSRs in analyzing the experimental data on 
7*(Q)7(g) vr'^-transition form factor. For ^ and <C m^, QCD factorization is valid 



only in the leading-twist approximation, so that the higher twists are important [231 ]. The reason is 
quite clear: if 0, one has to take into account the interaction of a real photon at long distances 
of the order of 0(1/y^?)- Then, in order to account for long-distance effects in perturbative QCD, 



one has to introduce the light-cone DA of the real photon. Instead of doing so, Khodjamirian 



4 



suggested to use the LCSR approach, which effectively accounts for the long-distances effects of 
the real photon using quark-hadron duality in the vector channel and a dispersion relation in 
q^. Schmedding and Yakovlev used this approach in analyzing the CLEO data on the 7*7 — > tt 
transition form factor at the level of the next-to-leading-order (NLO) accuracy of the perturbative 



QCD part of the LCSR 
We improved in 



21 



24 ] the NLO analysis of the CLEO data by taking into account the following 
key elements: (i) the NLO evolution for both Lp{x,Ql^^) and as(Qcxp) generalized to include 
heavy-quark thresholds more accurately; (ii) a relation between the "nonlocality" scale and the 
twist-4 magnitude ~ ^^g/^ was employed used to reestimate = 0.19 ±0.02 at A2 = 0.4 

GeV^; (iii) the possibility to extract constraints on {x~^)t^ from the CLEO data and compare them 
with those we derived before from NLC QCD SRs [15[ was exploited. 
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Fig. 3: Three 2cr- and Icr-contours (solid and dashed lines, correspondingly) of the admissible regions 
following from the analysis of the CLEO data for different values of A^: (a) = 0.4 GeV^; (b) A^ = 0.5 GeV^; 
(c) \l = 0.6 GeV^ 

n 

The results of our analysis in [21] are displayed in Fig. [3l Solid lines in all figures enclose the 
2(T-contours, whereas the Ic-contours are limited by dashed lines. The three slanted and shaded 

n 

rectangles represent the constraints on (02, 04) posed by the QCD SRs [151] for different values of 
Ag = 0.6, 0.5, 0.4 GeV^ (from the left to the right). Ah values are evaluated at ^? = 2.4 GeV^ 
after NLO evolution. 

We see that the CLEO data favor the value of the QCD nonlocality parameter = 0.4 GeV^. 
We also see from Fig. [3ljc) (and this conclusion was confirmed even with a 20% uncertainty of 
the twist-four magnitude — cf. Fig. Sl^a)) that the CZ DA (■) is excluded at least at the 4cj-level, 
whereas the asymptotic DA (♦) is off at the 3(T- level. At the same time, the BMS DA (X), and 
most of the bunch (the slanted green-shaded rectangle around the symbol IC), are inside the la- 
domain. The instanton-based Bochum (i^r) and Dubna (a) models are near but just outside the 



3o"-boundary and only the Krakow model 
the 2(T-boundary. 



251 ]. denoted in Fig. [D^a) by the symbol is close to 
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Fig. 4: (a): The results of the CLEO-data analysis for the pion DA parameters ((x )™''/3 — 1, evaluated at 
111 « 1 GeV^ (b): LCSR predictions for Q^i^^.^^^CQ^) for the CZ DA (upper dashed line), BMS- "bunch" 
(shaded strip), and the asymptotic DA (lower dashed line) in comparison with the CELLO (diamonds 2q ) 
and the CLEO (triangles 1221) experimental data, evaluated at /igy = 5.76 GeV^ with the twist-4 parameter 



value (5tw-4 



0.19 GeV^ 



2l| 



In Fig. IHa) we demonstrate the la-, 2a- and 3cj-contours (solid, dotted, and dashed contours 
around the best-fit point (+)), which have been obtained for values of the twist-4 scale parameter 

''Tw-4 



= [0.15 - 0.23] GeV^. As one sees from the blue dashed line within the hatched band 



(corresponding in this figure to the mean value of (2;~^)^^/3-l and its error bars) the nonlocal QCD 
sum-rules result with its error bars appears to be in good agreement with the CLEO-constraints on 
{x'^T^^^ at the lo--level. Moreover, the estimate is close to (x~^)™/3 - 1 = 0.24 ± 0.16, 

obtained in the data analysis of the pion electromagnetic form factor within the framework of 
a different LCSR method in 



27 



23] • These three independent estimates are in good agreement 
to each other, giving firm support that the CLEO-data processing, on the one hand, and the 
theoretical calculations, on the other, are mutually consistent. 
Another possibility, suggested in 



291 ] ■ to obtain constraints on the pion DA in the LCSR analysis 
of the CLEO data would be to use for the twist-4 contribution a renormalon-based model and to 
relate it to the parameters a2 and of the pion DA. Using this method, we obtained in [24 1 
renormalon-based constraints for the parameters a2 and 04 as shown in Fig. [6] in the form of a 
Icr-ellipse (dashed contour). 



III. DUET E791 DATA, PION FORM FACTOR AND CEBAF DATA 



Our findings are further confirmed by the E791 data [30(] on diffractive dijet vr -|- A-production. 
This is illustrated in Fig. [5ja). The main conclusion here is that all considered pion DAs are con- 
sistent with these data, with a slight preference for the BMS DA. Indeed, following the convolution 



procedure of |3l|, we found 



211 ] the following values of for the three types of pion DAs: 12.56 



6 



(asymptotic DA), 10.96 (BMS bunch), and 14.15 (CZ DA). 
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Fig. 5: (a): Comparison with the E791 data on the difFractive dijet production of the BMS "bunch" (shaded 
strip), the asymptotic DA (sohd line), and the CZ (dashed line) model, using the convolution approach of [3l|. 
(b): The scaled pion form factor calculated with the BMS bunch (shaded strip) and the asymptotic DA 



and 



15 



and the renormalization 



(dashed lines) including nonperturbative uncertainties from NLC QCD SRs 

scheme and scale ambiguity at the O(cks) level |33|. The experimental data are taken from [Sj] (diamonds) 



351] (triangles). 



It is worth mentioning the results of our analysis of the pion electromagnetic form factor using 
the NLC-based pion DA and analytic perturbative QCD [33] . These results are in excellent agree- 
ment with the CEBAF data on the pion form factor, as illustrated in Fig. [5^b), where the green 
strip includes both the NLC QCD SR uncertainties, generated by the bunch of the allowed pion 
DAs, and by the scale-setting ambiguities at the NLO level of QCD perturbation theory. 

From the phenomenological point of view, the most interesting result here is that the BMS 
pion DA [15] (out of the "bunch" of similar doubly-peaked endpoint-suppressed pion DAs) yields 
predictions for the electromagnetic form factor of the pion, which are very close to those obtained 
with the asymptotic pion DA. Conversely, we see from Fig. [5][b) that a small deviation of the 
prediction for the pion form factor from that obtained with the asymptotic pion DA (dashed lines) 
does not necessarily imply that the underlying pion DA has to be close to the asymptotic profile. 
Much more important is the behavior of the pion DA in the endpoint region x 0,1. For more 



details, we refer to 33|] 



IV. NEW LATTICE DATA AND PION DA 



Rather recently, new high-precision lattice measurements of the second moment of the pion DA 
= JoC^^ ~ ^)'^^Tr{x) dx appeared [s^]- Both cited groups extracted from their respective 
simulations values of 02 at the Schmedding-Yakovlev scale /igy around 0.24, but with different 
error bars. 
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It is remarkable that these lattice results are in striking agreement with the estimates of oo from 



both the NLC QCD SRs |15] and also from the CLEO-data analyses — based on LCSR- 
as illustrated in Fig. where the lattice results of [32] are shown in the form of a vertical strip, 
containing the central value with associated errors. Remarkably, the value of 02 of the displayed 
lattice measurements (middle line of the strip) is very close to the CLEO best fit in 21[ (+). 



V. IMPROVED MODEL FOR NLCS AND PION DA 

The quark-gluon-antiquark condensates are usually parameterized in the following form 

(^(0)7m gW) = {y^,x,-g^,{yx))M,{x^y^,{y-xf) 

+ iVfiVu - gtiuy^)M2ix'^,y^,{y - xf) , 
(9(0)757m = ie^,y,M3{x^,y\{y-xf), 

00 

where Mi{x\y\z'^) = Ai [[[dadp d-f fi{a, e("'''+'^y'+^'')/'^ with ^,2,3 = (-|,2, |). The 





minimal model of the nonlocal QCD vacuum suggests the following Ansatz 

(a, /3, 7) = 5 (a - A) <5 (/? - A) ^ (7 - A) (9) 

with A = and faces problems with the QCD equations of motion and the gauge invariance of 
the 2-point correlator of vector currents. In order to fulfil the QCD equations of motion exactly 
and minimize the non-transversity of V — V correlator, an improved version of the QCD vacuum 
model was suggested in [3j] with 



/r"" (a, 13, 7) = (1 + X,d, + Yidy + Y,d,) 5 (a - xA) 6 {(3 - yA) 6 (7 - zA) , (10) 



where A = ^A^ and 



Xi = +0.082; ^2 = -1.298; ^3 = +1.775; x = 0.788; (11a) 
Yi = -2.243; ^2 = -0.239; = -3.166; y = z = 0.212. (lib) 

Then, the NLC sum rules ([5]) gives rise to a modified "bunch" of two-parameter pion DA models 
([6]) at /i^ = 1.35 GeV^ 32] • The coordinates of the central point ■♦■ are a2 = 0.268 and 04 = —0.186. 
These values correspond to (x~^)^"°'^^ = 3.24+0.20, which is in agreement with the result provided 
by an independent sum rule, viz., {x~ 

-i)SR = 3.40 + 0.34. The allowed values of both bunch 
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Fig. 6: The results of the CLEO-data analysis for the pion DA parameters 02 and 04, evaluated at /igy = 
5.76 GeV^, are shown. The lattice results of [37| are displayed as a shaded area, whereas the renormalon- 
based lu-ellipse of 2^ is denoted by the green dashed line. The small slanted rectangle on panel (a) 
represents the corresponding results for the BMS- "bunch" , whereas the modified, somewhat larger, slanted 
rectangle corresponding to the improved Gaussian model of QCD vacuum, is shown on panel (b). 

parameters 02 and 04 after NLO-evolution to = 5.76 GeV^ are shown in Fig. [6] in the form of 
shaded slanted rectangles around the central points K and 

We emphasize in this context that the BMS model [l^, shown in Fig. [6)^a) by the symbol K, is 
inside the allowed region extracted from the improved QCD vacuum model. This means that all 
the characteristic features of the original BMS bunch are also valid for the improved bunch. Again 
the NLC-based model DAs are end-point suppressed, though they are doubly humped. 

We see from Fig. El^b) that the improved bunch 
with the recent lattice results 



32( 1 appears to be in somewhat better agreement 



371], shown in the form of a vertical strip, containing the central 
value with its associated errors. Remarkably, the value of 02 of the displayed lattice measurements 
(middle line of the strip) is very close to the central point of the improved bunch (>■), whereas the 
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Fig. 7: Evolution of (^^}^(^^) (left panel) and (C^)Tr(Ai^) (right panel) with G 1 — lOGeV^. The green strip 
on both panels corresponds to the unified results of the QGD Sum Rules with NLCs (minimal and improved 
models of the QCD vacuum) . On the left panel, we also show the lattice results with their corresponding 



error-bars: 



361, ♦ [371, ■ 129|- 
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whole admissible region (slanted rectangle) 84] is inside the strip^ and also inside the standard 
CLEO la-ellipse. Furthermore we see from this figure that the higher the precision of the lattice 
simulations, the closer they are to the results for (C^)7r(/U^) of the NLC-based QCD SRs. It remains 



to be seen whether this agreement wil 
will become possible on the lattice [4C|. 



persist also for the {C'^)n{fJ''^) moment, once it is calculation 



VI. PION DA AND DRELL YAN ttN PROCESS 

The DY process is the dominant mechanism to produce lepton pairs with a large invariant 
mass in hadronic collisions, like vr^A^ scattering. In this model a massive muon pair is created 
through the electromagnetic annihilation of an antiquark from the beam pion and a quark from 
the nucleon target, as depicted on the left panel of Fig. El 




Fig. 8: Left panel: Graphical representation of the Drell-Yan process -k^ N fi^fji^X. Right panel: 
Angular definitions of the Drell-Yan process in the center of mass frame of the produced massive lepton 
pair. The axis zgj denotes the pion direction in the Gottfried-Jackson (GJ) frame. 

Here pu = xuP is the momentum of the annihilating antiquark from the pion. Typical values 
of the kinematical parameters s = {p-,^ + Pn)'^, = <f, and = —q^, see Fig. (H in the case of 
the FNAL experiment E615 are: s = 500 GeV^, = - 70 GeV^ and p = Q^/Q = - 0.5. As 

— > 1, becomes large and far spacelike and, therefore, it is sufficient to consider the u-quark 
to be nearly free and on-shell: = (no transverse momenta). On the right panel of Fig. 
[8] we show the angular-distribution parameters 9, the polar angle measuring the direction in 
the Gottfried-Jackson system of axes, and <j), the azimuthal angle between the 7r~/u and tt^N 
planes in the lepton rest frame. 



^ This statement is valid for the BMS "bunch" as well, as can be seen from Fig. ^a). 
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Fig. 9: Results for the angular distribution parameter ^ as a function of Xu = for different values 
of p = Qt/Q- The green strip corresponds to the BMS bunch of the pion DAs [l^ with the solid line 
representing the BMS DA, while the dotted solid line shows the result for the asymptotic DA, and the 
dashed line denotes the prediction for the endpoint-dominated CZ DA. The data are taken from 

For the DY reaction with an unpolarized target, the angular distribution of the in the pair 
rest frame can be written in terms of the kinematic variables A, /i, i' as follows 



dQ'^dQ'i.dx ^ d cos 



oc iV(x, p){l + X cos^ 61 + /i sin 26 cos + ^ sin^ 6 cos 2(/>^ . (12) 



Adopting this convolution procedure, we found 



431] the results presented in Fig. [D We see that 



the agreement of the chosen pion DA model with the unpolarized E615 (FNAL) data depends on 
the value of the parameter p. It seems that these data cannot make a clear distinction in favor of 
one particular pion DA. On the other hand, for the asymmetry of the polarized DY tt^ N process 
we found (using the convolution procedure of [44] ) the results displayed in Fig. [TUl We may come 
to the conclusion that the asymmetry A{(j), xl) can be used to discriminate different proposed pion 



DA models, provided the value of p can be fixed by experiment 



43] 




. 4 

X]^ 0.6 



Fig. 10: 3D-plots of the azimuthal asymmetry A{4',XL,p = 0.3) at = 16 GeV^ and s — 100 GeV^ for 
three different choices of the pion DA: Asymptotic (left), BMS model (center), and CZ model (right). 
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VII. CONCLUSIONS 



Let us conclude with the following observations: 
(i) The NLC QCD SR method for the pion DA gives us admissible bunches of pion DAs for each 
value of Xg > 0.4 GeV^. (ii) The NLO LCSR method produces new constraints on the pion DA 
parameters (02 and 04) in conjunction with the CLEO data, (iii) Comparing the results of the 
NLC SRs with the new CLEO-data constraints, allows us to fix the value of the QCD vacuum 
nonlocality to be Xg ~ 0.4 GeV^. (iv) The bunch of pion DAs from NLC QCD SRs agrees well 
with the E791-data on the diffractive dijet vr + A-production, with the JLab F(pi) data on the pion 
electromagnetic form factor, and with the latest high-precision lattice data, (v) Taking into account 
the QCD equations of motion for the NLCs and the transversity of the vacuum polarization, allows 
us to shift the pion DA bunch just inside the la-ellipse of the CLEO-data constraints, (vi) In the 
polarized Drell-Yan 7r~A^ process one might be able to discriminate among different pion DAs, 
once the value of p will be known by measurements. 
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